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ABSTRACT 

We observed two faint tidal dwarf galaxies (TDGs), NGC 5291 N and NGC 5291 S with the Infrared 
Spectrograph on the Spitzer Space Telescope. We detect strong polycyclic aromatic hydrocarbon 
(PAH) emission at 6.2, 7.7, 8.6, 11.3, 12.6 and 16.5 /im, which match models of groups of ^ 100 
carbon atoms with an equal mixture of neutral and ionized PAHs. The TDGs have a dominant warm 
~ 140 K dust component in marked contrast to the cooler, 40 — 60 K dust found in starburst galaxies. 
For the first time we detect the low-J rotational lines from molecular hydrogen. Adopting LTE there 
is ~ 10^ M0 of ~ 400 K gas, which is < 0.1 % of the cold gas mass inferred from ^^CO (1-0) 
measurements. The combination of one-third solar metallicity with a recent, < 5 million year, epsiode 
of star formation is reflected in the S and Ne ratios. The excitation is higher than typical values for 
starburst galaxies and similar to that found in BCDs. Using the Infared Array Camera we identify an 
additional 13 PAH-rich candidate TDGs. These sources occupy a distinct region of IRAC color space 
with [3.6] — [4.5] < 0.4 and [4.5] — [8.0] > 3.2. Their disturbed morphologies suggest past merger 
events between companions, for example, NGC 5291 S has a projected 11 kpc tail. NGC 5291 N and 
S have stellar masses of (1.5 and 3.0) x 10* Mq, which is comparable to BCDs, though still roughly 
10% of the LMC's stellar mass. The candidate TDGs are an order of magnitude less massive. This 
system appears to be a remarkable TDG nursery. 

Subject headings: galaxies: dwarf — galaxies: interactions — individual ( NGC 5291, NGC 5291 S, 
NGC 5291 N) — galaxies: formation — infrared: galaxies 
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1. INTRODUCTION 

In addition to triggering starbursts and active galactic 
nuclei, mergers of dusty, gas rich disk galaxies frequently 
lead to the formation of tidal tails stretching >100 kpc 
from the site of the collision (Toomre & Toomre 1972; 
Schweizer 1978; Sanders & Mirabel 1996). These struc- 
tures tend to be HI rich with blue optical colors, reflect- 
ing both their origin in the outer spiral disk and on-going 
star formation (van der Hulst 1979; Schombert Wallin & 
Struck-MarceU 1990, Mirabel, Lutz & Mazza 1991; Hib- 
bard & van Gorkom 1996). Zwicky (1956) proposed that 
dwarf galaxies might form out of self-gravitating clumps 
within tidal tails, and indeed, concentrations of gas and 
star forming regions, with dwarf galaxy size HI masses 
and optical luminosities, are commonly found. Computer 
simulations of merging galaxies typically give timescales 
of ^ 10* years for the formation of condensations at the 
tips of the ^ 100 kpc tails, indicating that any young 
stellar population has been formed in-situ and has not 
simply been stripped from the host galaxy disk. 

The cataloging of Tidal Dwarf Galaxies (TDGs) is still 
in its infancy. Kinematic observations alleviate potential 
contamination from foreground and background objects. 
However projection effects can give the appearance of 
a bound object to unbound condensations in the debris 
left over from galaxy interactions. A bona flde TDG 
should be self-gravitating, rotating, long-lived and tidal 
in origin. Tidal origin and longevity can only be inferred 
from comparison of the source morphology and kinemat- 
ics with simulations. It is often impossible to fully disen- 
tangle tidal effects from other dynamical effects such as 
ram pressure stripping. Whether or not most candidate 
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TDGs are subject to tidal disruption and re-absorption 
by the merger remnant on ^Gyr timescales is still an 
open question (e.g., Hibbard & Mihos 1995; Elmegreen, 
Kaufman & Thomasson 1993). The formation of both 
TDGs and transient unbound star-forming regions will 
enrich the outer regions of the inter-galaxy medium. 

The study of dwarf irregular proto/young galaxies 
formed either via tidal interactions between the parent 
galaxies or from ram-sweeping of debris material, and 
collectively called TDGs, may be an important part of 
galaxy formation and evolution. TDGs located at the 
tips of optically faint but HI rich tails represent new star 
forming systems largely free of pre-existing stellar popu- 
lations, and represent useful templates of star formation 
in the early universe, albeit with higher metallicity (Due 
& Mirabel 1998, hereafter, DM98). More importantly, 
dwarf galaxies are by far the most common galaxy type 
in the current epoch, accounting for 38 out of 41 galaxies 
in the Local Group alone (Mateo 1998). Work by Huns- 
berger et al. (1996) suggests that TDGs may make up 
at least 1/3 and possibly 1/2 of the dwarf galaxy pop- 
ulation in compact groups i.e., a potentially significant 
galaxy formation mechanism in the current epoch. One 
problem with this scenario is dark matter. Like mas- 
sive spirals, dwarf galaxies appear to be dominated by 
dark matter. If dark matter is non-baryonic it is hard to 
form a halo around a dwarf galaxy that has been formed 
via tidal interaction. If the dark matter is baryonic i.e., 
cold molecular gas (Pfenniger, Combes & Martinet 1994) 
then the TDG/dwarf association remains a valid possi- 
bility. Both the nature of dark matter and whether the 
dynamical masses of TDGs infers its presence are still 
open questions. 

Studies of TDG properties may help us better under- 
stand the dwarf galaxy population as a whole, such as 
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what fraction are truly primordial building blocks of mas- 
sive galaxies left over from the epoch of galaxy forma- 
tion, and what fraction may be TDGs from high and 
intermediate red-shifts. This is especially relevant as the 
merger rate was likely much higher at z > 1 than it is 
today. To date the TDG sample is still too small to 
statistically assess the importance of this population of 
objects. That assessment will be made in due course. 
Whether rare or ubiquitous, these objects warrant de- 
tailed study. Here we utilize the unique capabilities of 
the Spitzer Space Telescope (Werner et al. 2004a) to 
study the NGC 5291 system, which is known for its re- 
markably rich proto/young galaxy nursery. 

NGC 5291 is a disturbed elliptical/lenticular galaxy 
located at the western edge of the cluster Abell 3574. 
Pedersen et al. (1978) briefly discussed deep images of 
the system taken with the ESO 3.G m telescope. They 
highlighted the presence of many optical knots, which 
they suggested were extragalactic HII regions. Long- 
more et al. 1979 (hereafter, L79) published a much 
more detailed investigation. Using deep Schmidt plates 
with filters approximating B (AA 3950 - 5400 A Schott 
GG395 mter/III-aJ emulsion) and U (AA 3200 - 3900 
A Schott UGl filter /Ill-a J emulsion) passbands they 
found 24 blue knots extending across an arc centered on 
NGC 5291 and obtained spectra of the eight brightest. 
They identified the knots as giant extragalactic HII re- 
gions, but noted that both the projected linear size and 
optical luminosity of the largest knots were greater than 
the Large Magellanic Cloud. Coincidentally, Maza et 
al. (1991) also obtained spectra of the brightest opti- 
cal knot, which was included in an object prism survey 
of HII galaxies. Using the Parkes radio telescope, L79 
also discovered a large HI mass (~10^^ M©) associated 
with NGC 5291 and its blue knots though slightly dis- 
placed to the west. Malphrus et al. (1997, hereafter, 
M97) mapped the HI 21 cm emission using the VLA. 
The HI morphology resembles a large fragmented ring 
with a continuous eastern arc stretching ~ 4 arcmin (72 
kpc) on either side of NGC 5291. The eastern half of the 
ring contains most of the neutral hydrogen as well as the 
highest HI surface densities, and is associated with both 
NGC 5291 and its knots of star formation. The origin of 
this giant HI ring is still unclear. NGC 5291 has a dis- 
turbed companion galaxy called the "Seashell" , however 
it is unlikely that this galaxy is linked to the HI ring's 
formation as it is both gas poor and has a much higher 
radial velocity (~600 km s~^) relative to NGC 5291, indi- 
cating a high fly-by velocity (L79, > 400 km s-\ DM98). 
L79 discussed several models for the origin of NGC 5291. 
They noted that the optical morphology is similar to the 
Antennae (NGC 4038/39) and that there are two other 
galaxies (# 435 and 446 in Richter 1984) within 150 kpc 
projected distance from NGC 5291 that may be respon- 
sible for a tidal interaction. However M97 rule out the 
interplay for both galaxies. They detected Richter 435, 
which hes ~ 4' southwest of NGC 5291, in HI, but with 
a heliocentric velocity ~ 250 kms^^ lower than the HI 
ring. M97 found that, like the Seashell, there is little HI 
associated with the second candidate galaxy Richter 446, 
~ 5' northeast of NGC 5291. However, L79's preferred 
model for the system is that a pre-existing giant HI disc 
has undergone compression and possible ram pressure 



sweeping by the intra-cluster medium as the system falls 
into the cluster from the far side. L79 suggested that the 
knots could eventually form a population of blue dwarf 
irregular galaxies. The combined optical and VLA ob- 
servations led M97 to revive the tidal model and they 
conclude that the system is probably both interacting 
and likely being ram-swept as it moves through the intra- 
cluster medium. 

However mysterious its origin the NGC 5291 system is 
a remarkable galactic nursery. M97 investigated whether 
the knots are self-gravitating and will remain bound 
against the tidal force of the parent galaxy to become 
dwarf irregular galaxies. They compared the HI mass 
for each knot to the binding mass. They argued that for 
the two most massive HI clumps, NGC 5291 N (hereafter, 
TDG-N), and NGC 5291 S (hereafter, TDG-S) the com- 
bined mass of stars and gas will result in self-gravitating 
knots. M97 also calculated the radius at which material 
in a knot would be disrupted by the tidal force with the 
parent galaxy and compared it to the projected radius 
of the HI knots. The majority of the knots, including 
TDG-N and TDG-S, appear to be tidally stable. Us- 
ing Ha Fabry-Perot observations Borunaud et al. (2004) 
were able to determine the rotation curve for TDG-N, 
clearly showing that it is self-gravitating and rotating. 
They concluded that higher resolution observations are 
required to assess the presence of dark matter. 

Before the launch of the Spitzer Space Telescope, mid- 
infrared spectroscopic studies of objects as faint as TDGs 
were not feasible. We have taken advantage of the un- 
precedented sensitivity of the Spitzer Space Telescope to 
obtain IRS spectra of the two optically brightest TDGs 
in the NGC 5291 system, TDG-N and TDG-S at the 
northern and southern-most tips of the high surface den- 
sity HI arc. The mid-infrared region of the spectrum 
is rich in atomic and molecular emission features from 
the star forming ISM, including fine structure lines like 
[Ne II], polycyclic aromatic hydrocarbons (PAHs) and 
warm molecular hydrogen. This will allow comparisons 
with similar mid-infrared studies of other star forming 
galaxies from dwarf irregulars to massive spirals. 

We have also imaged the eastern half of the HI ring, in- 
cluding the associated proto/young irregular galaxy pop- 
ulation, NGC 5291 and the "Seashell", with the Infrared 
Array Camera (Fazio et al. 2004, IRAC). We assess the 
relative contributions from both cool (evolved and low 
mass main sequence) stars and from star formation as 
traced by PAH emission in the 8.0 /im band. The result- 
ing mid-infrared spectral energy distribution (SED) and 
colors will allow useful comparisons with other galaxy 
types. 

In §2 we detail the observations and data reduction. 
§3 begins with the results and discussion of the the mid- 
infrared spectral properties of TDG-N & TDG-S. We 
compare the stellar and non-stellar (i.e., PAH) morphol- 
ogy of NGC 5291, the "Seashell", and the knots next 
in §§3.2 — 3.4. In §3.5 we select a population of can- 
didate TDGs by matching the IRAC SEDs to galaxy 
templates. We investigate the mid-infrared color-color 
and color-magnitude relations in §3.6. We estimate the 
TDG stellar masses in §3.7 and present our conclusions 
in §4. We adopt a luminosity distance of 62 Mpc to the 
NGC 5291 system. This distance is calculated using the 
red-shift in NED with H^ = 71 kms'^Mpc'S Om = 
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0.27, r^A = 0.73. & Ofc = 0, and does not reflect any 
uncertainty in the distance due to proper motion within 
the group. This gives a hnear scale of 1" = 0.3 kpc. 

2. OBSERVATIONS & DATA REDUCTION 
2.1. IRS Spectroscopy 

IRS observations were taken of both TDG-N and TDG- 
S. The IRS low resolution data were obtained with 
Short-low (IRS-SL), which operates between 5.2 and 7.7 
/im (IRS-SL2) and 7.4 - 14.5 ^iul (IRS-SLl). IRS-SL has 
a resolving power of 64 < < 128. The high resolution 
data were obtained using Short-high (IRS-SH), which en- 
compasses the range 9.9 - 19.6 /zm , and Long-high (IRS- 
LH), which spans 18.7 - 37.2 fj,m. The two high resolution 
spectrometers have a resolving power, ~ 600. 

Observations were made in the IRS Staring Mode AOR 
with a high accuracy blue peak-up using a star from the 
2MASS catalogue (Cutri et al. 2003). The staring mode 
AOR integrated for half the time listed in Table 1 at each 
of the two nominal nod positions on each slit. The spec- 
tral data were processed as far as the un-flat-fielded two 
dimensional image using the standard IRS Sll pipeline. 
See chapter 7 of the Spitzer Observers Manual (SOM)^ 
for further observing mode and pipeline details. 

Both TDG-N and TDG-S have extended emission at 8 
^m with respect to the IRS slits (3.6" IRS-SL, 4.7" IRS- 
SH and 11.1" IRS-LH). The spectra were extracted and 
sky subtracted using the extended source routines in the 
SMART analysis package (Higdon et al. 2004). The 
IRS-SH and IRS-LH date were extracted using the full 
aperture and the IRS-SL spectra were extracted using a 
column of 20" fixed width. The spectra were flat-fielded 
and flux-calibrated by extracting and sky subtracting un- 
flat-fielded observations of the calibration stars HR 6348 
(IRS-SL) and ksi Dra (IRS-SH, IRS-LH), and dividing 
these data by the corresponding template (Cohen 2003) 
to generate a 1-dimensional relative spectral response 
function (RSRF). The RSRF was then applied to the 
TARGET observations to produce the final spectra. The 
residual sky was subtracted from the IRS-SL data using 
the off-source observations, which are part of the Star- 
ing mode AOR. A Zodiacal light model from the Spitzer 
Planning Observations Tool (SPOT) was scaled to the 
IRS slits and subtracted from the IRS-SH and IRS-LH 
data. An aperture correction of 0.85 was applied to the 
IRS-SL spectra as the sources are extended. In addition, 
the IRS-LH TDG-N and TDG-S spectra were scaled by 
0.30 and 0.26, respectively in order to stitch to the IRS- 
SH spectra. The scale factors are close to the simple 
geometric aperture correction for the differing IRS-SH 
and IRS-LH slits which is 0.215. However, these sources 
are not uniformly extended and the two apertures are 
orientated ^ 85° with respect to each other. 

2.2. IRAC Imaging 

IRAC data in all four bands were acquired using 
a two point map centered on NGC 5291 (q!J2ooo = 
13'*47™23.00^ (5j2ooo = -30°25'30.0"). The array was 
aligned to celestial north and stepped by 180". Three 
3 s frames were taken in each channel using a 12 posi- 
tion Reuleaux dither pattern with a medium scale fac- 



: //ssc . spitzer . caltech. edu/documents/som/ 



tor, resulting in a total integration time of 432 s. Due to 
the footprint of the IRAC arrays the resulting map has 
a total area of 110.8 arcmin^ with a sub-region of 59.1 
arcmin^ sampled by all 4 channels. In this paper we will 
only discuss the properties of sources located in this sub- 
region. The IRAC data were processed through version 
Sll. 0.2 of the SSC-IRAC pipehne. For more observing 
mode and pipeline information see chapter 6 of the SOM. 

The IRAC was designed to map both interstellar emis- 
sion and starlight. The 3.6 & 4.5 /im bands primarily 
trace the stellar mass distribution, as the emission at 
these wavelengths is well approximated by the Rayleigh- 
Jeans limit of blackbody emission for both early- and 
late-type stars. At 8 /j,m the mid-IR traces dust emis- 
sion, in particular the 7.7 & 8.6 /im PAH features from 
the interstellar medium. It is straight forward to subtract 
the stellar emission from the 8 /im image to produce the 
PAH/dust image shown in Figure 1. The method is out- 
lined in Pahre et al. (2004). The 3.6 and 4.5 /im images 
are scaled to match the theoretical colors of MO HI stars 
([3.6] - [4.5] = -0.15, [4.5] - [5.8] = +0.11 and [5.8] - [8.0] 
= +0.04) at 8 /xm, and subtracted from the 8 /im image. 

SExtractor (Bertin & Arnouts 1996) was used to gen- 
erate a catalog of sources for each IRAC channel. We 
then selected all sources that were identified in all four 
channels with a positional uncertainty of < 2 ". SEx- 
tractor was also run on the dust image to find an ad- 
ditional 14 sources which were not detected in all four 
channels. For these sources SExtractor was re-run on 
the 3.6 through 8.0 /im images using the source posi- 
tions determined from the 8 /im image. Table 4 lists the 
flux densities in an 8 pixel (9.8 ") diameter aperture. A 
16 pixel (19.5 ") diameter aperture was used for a few 
sources, including TDG-N and TDG-S. The flux densi- 
ties were corrected to a 20 pixel radius aperture using 
the following scale factors 1.086 & 1.017 (3.6 /im), 1.098 
& 1.018 (4.5 /xm), 1.097 & 1.020 (5.8 /xm) and 1.118 & 
1.026 (8.0 /im) for the 8 and 16 pixel apertures respec- 
tively, as determined by the IRAC team. The auto-flux 
returned from SExtractor is included for objects that are 
extended beyond the 19.5 " circular aperture, such as 
NGC 5291 and the Seashell. The Vega zero magnitude 
flux densities in the four IRAC bands are 277.5, 179.5, 
116.6 and 63.1 Jy at 3.6, 4.5, 5.8 and 8.0 /im respectively. 
Stellar sources were removed from the catalog if they 
matched one of the following criteria: 2MASS (J - K) 
< 0.5 (Eisenhardt 2004), identified as stellar sources in 
the Hubble Guide Star Catalogue Version 2.2, or through 
visual inspection of the IRAC images. 

3. RESULTS AND DISCUSSION 

3.1. Mid-Infrared Spectral Properties of TDG-N & 
TDG-S 

The top of Figures 2 & 3 show the IRS-SL spectra of 
TDG-N and TDG-S. The broad emission features from 
PAHs at 6.2, 7.7, 8.6 and 11.3 12.6 /im are clearly present. 
Bright fine structure lines of [S iv] 10.51 /im and [Neil] 
12.81 /im emission lines are also apparent in this low res- 
olution spectrum, as is a weak rotational 0-0 S(3) 9.66 
/tm line from molecular hydrogen. Figure 4 shows the 
IRS-SL spectrum of TDG-N overlaid with the IRS-SL 
spectrum of the reflection nebulae NGC 7023 (Werner 
et al. 2004b) and the ISO-SWS spectrum of M 82's nu- 
cleus (Lutz et al. 1998), smoothed to the IRS-SL reso- 
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lution. The overlay spectra are normalized to TDG-N's 
flux density at 7 /im. Here we arc using NGC 7023 to 
characterize a low ionization star forming region as it is 
illuminated by HD 200775, a B2Ve star. M 82 is used to 
represent the mid- infrared spectrum of a starburst. The 
three spectra are remarkably similar in shape and only a 
detailed analysis will reveal their differing environments. 
To assess whether the PAHs in TDG-N and TDG-S are 
similar to those in M 82 we fit the continuum emission in 
the IRS-SL spectrum with the thermal radiation from a 
distribution of graphite and silicate grains and the PAHs 
with a series of Drude profiles. Table 2 lists the resulting 
PAH fluxes and equivalent widths. The relative strength 
of the 11.3/7.7 fim PAHs in TGN-N (TDG-S) is 0.22 
t 0.03 (0.23 t 0.05) and 6.2/7.7 ^im PAHs 0.35 t 0.06 
(0.33 1 0.11). These values are consistent with groups of 
~ 100 carbon atoms in an approximately equal mixture 
of neutral and ionized PAHs (see Figure 16 in Draine 
& Li 2001). For comparison, on the same Figure, the 
PAH ion fraction is close to unity in M 82 with similar 
sized clumps of C atoms. In addition to the PAH fea- 
tures observed in the IRS-SL spectra, we also detect PAH 
emission at 13.5 and 16.5 /um in the IRS-SH spectra, the 
16.5/7.7 /im PAH ratio is 0.02^0.01 and 0.10t0.05 in 
TDG-N and TDG-S. 

Figures 2, bottom & 3, bottom show the IRS-SH & 
IRS-LH spectra. The fine structure neon and sulfur lines 
are prominent as well as weak 0-0 S(l) and 0-0 S(2) emis- 
sion from H2. The spectra flatten out around A >20 /xm. 
A similar flattening is seen in the extremely low metallic- 
ity blue compact dwarf galaxy SBS0335— 052 (Houck et 
al. 2004b). Using SMART we fit a diluted black-body to 
the continuum over the range 15 — 35 fim . Using a solid 
angle equal to the IRS-SH aperture of 1.35 x 10~^° sr 
gives a temperature T = 139 1 6 K. The derived optical 
depth at 0.55 /zm (ro.ss^m) is effectively zero, and a = 
1.2 't 0.2 where, r = To.55^m(0.55/A)". A similar result 
is obtained for TDG-S with a temperature T = 142 1 31 
K. A dominant cooler 40 — 60 K dust component that is 
present in most starbursts and spiral galaxies is absent. 

Figures 5 & 6 show the line profiles from the IRS obser- 
vations. The line fluxes and equivalent widths for both 
TDGs arc given in Table 3. No correction has been made 
for reddening as the optically derived extinction is low, 
typically, Ab < 1 mag. (L79, Pena, Ruiz & Maza 1991, 
DM98). Geometric effects can also make any extrapola- 
tion from the optical to the infrared highly uncertain. In 
the following paragraphs we use the line strengths and 
ratios to derive physical properties of the ionized and 
molecular ISM in the TDGs. 

Both the Ne and the S line ratios can be used to con- 
strain the properties of the underlying starburst. The 
starburst sample of Thornley et al. (2000) have typi- 
cal [Ne iii]/[Ne 11] ratios between 0.05 and 1.0. For ex- 
ample, M 82 has a neon ratio of 0.18. The exceptions 
to this range are the two low metallicity dwarf galaxies, 
NGC 5253 and II Zw 40, which both have one-fifth so- 
lar metallicity. These have neon ratios of 3.5 and 12, 
respectively. Lowering the metallicity produces "hotter" 
main sequence stars for a given mass, and the radiation 
is harder due to reduced line blanketing and blocking, 
hence the neon ratio increases. The TDG-N and TDG-S 
have neon ratios of 2.4 lO.l and 1.4 1 0.1 respectively, 



which is consistent with their one-third solar metallicity 
(L79, Pena, Ruiz & Maza 1991, DM98). Figure 10 from 
Thornley et al. (2000) gives an age for the most recent 
episode of star formation of ~ 5 million years for these 
neon ratios for a galaxy with one-fifth solar metallicity 
and log U = — 2.3, where U is the number of ionizing 
photons at the surface of the nebula per hydrogen atom. 
This result is consistent with the starburst age derived 
from optical spectra by DM98. Both TDG-N and TDG- 
S have slightly higher metallicity, as well as fewer ioniz- 
ing photons (i.e., log U ~ —3), as inferred from model 
excitation diagrams in Pena, Ruiz & Maza (1991, for 
TDG-S we used emission lines from DM98). Both these 
effects will result in a younger burst age using the figure 
from Thornley et al. (2000). This is consistent with the 
prediction that the stars have formed in-situ and were 
not simply dragged from the outer regions of the parent 
galaxy during a tidal interaction. 

The [Siv]/[Siii] line ratio is another excitation diag- 
nostic. In M 82 this ratio is 0.05 (Verma et al. 2003), 
compared to 0.48 t 0.02 and 0.23 t 0.05 in TDG-N 
and TDG-S, respectively. Verma et al. (2003) plot- 
ted the log of the [Neiii]/[Neii] against the log of 
the [Siv]/[Siii] line ratio for a sample of twelve star- 
burst and blue compact dwarf (BCD) galaxies. Star- 
bursts tend to occupy the bottom left quadrant of the 
plot, having low excitation. The upper right quadrant 
is filled with BCDs and galaxies exhibiting Wolf-Rayet 
features. The corresponding ratios for TDG-N and TDG- 
S place them in the upper right quadrant with moder- 
ate excitation. The [Slv]/[Slll] line ratio in TDG-N 
is close to the value for the overlap region in the in- 
teracting galaxy pair NGC 4038/NGC 4039, but the 
[Ne III] /[Nell] ratio is a factor two higher (see Figure 
4 in Verma et al. 2003). The ratios for TDG-S are lower. 
The [Siv]/[Siii] line ratio is similar to the starburst 
galaxy, NGC 7714 [Siv]/[Siii] = 0.20, but again the 
[Ne lll]/[Ne 11] ratio is a factor two higher (Brandl et al. 
2004). 

The [Slll] 33/18.7 fim ratio is a sensitive function of 
electron density in the range 300 - 30,000 cm~^. The 
[Siii] 33 /um line is detected in TDG-N though with 
low signal-to-noise. The [S III] ratio is 0.7 ^ 0.3, which 
implies an electron density of 600 < < 3000 cm~^ 
(Figure 1 in Rubin et al. 1994). L79 derived a limit 
of < 10'^ cm~^ from their optical spectrum, which is 
consistent with our result derived from infrared lines. 

We detect weak H2 emission in both TDG-N & TDG-S. 
The 0-0 S(2) 12.28 iim and 0-0 S(l) 17.03 /xm rotational 
emission lines have low signal-to-noise (SNR > 3) in the 
IRS-SH spectra, but the 0-0 S(3) 9.66 /xm is also visible 
in the IRS-SL spectra (sec top of Figures 2 & 3), and 
all lines, except for the 0-0 (S2) line in TDG-S, are de- 
tected in the individual spectra from both nod positions, 
adding weight to the detection of H2 . The average line 
profiles are shown in Figures 5 & 6. To derive the mass of 
warm molecular hydrogen we assume that the emission 
is optically thin. The critical densities of the low J levels 
are relatively low (ucr < 10^ cm~^) and we assume that 
the populations are in LTE. Adopting an ortho to para 
ratio of 3 we show the excitation diagrams in Figures 7. 
These are simply the natural logarithm of the number of 
molecules divided by the statistical weight in the upper 
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level of each transition versus the energy level. For a 
single component model the data lies on a straight line 
and the excitation temperature is the reciprocal of the 
fit to the slope. The mass is derived from the S (1) line 
luminosity and the excitation temperature. In the fol- 
lowing equations we outline the derivation for the mass 
of molecular hydrogen. The energy of a given level, J, is 

Ej = 85kJ{J+l) (1) 
where k is the Boltzmann constant. The total mass is 



MTotal — g^'-^o- 

Here Mq is the mass of gas in the ortho state and 



(2) 



Mo = mH^Nr (3) 

where mjjj is the molecular mass of H2 and Nt is the 
total number of molecules. 



JJ 

Nj is the number of molecules for the J*'' state 



Nj = 



m 

AjAEj 



(4) 



(5) 



L(J) is the line luminosity, Aj is the Einstein A coeffi- 
cient, AEj = hiyj, where h is Planck's constant and vj 
is the frequency of the emission line, f j is the partition 
function for the J*'' state 



fj = 



gje 



"^J'. ortho f}J' 



(6) 



where T is the excitation temperature and gj is the 
statistical weight for a given state. 

For TDG-N and TDG-S the temperature is 460 t 75 K 
and 400 i 40 K. The value is slightly higher than the av- 
erage vahie of ^ 350 K measured by us for a large sample 
of ULIRGs. Assuming the primordial helium mass frac- 
tion is ~ 0.25 (e.g. Olive & Skillman 2004), we derive 
a total mass of warm gas of (1.0 0.7) and (1.5 0.6) 
xlO^ Mq in TDG-N and TDG-S, respectively. The 0-0 
S(0) limit is consistent with this temperature, but does 
not rule out a colder, say 155 K component, which would 
imply a mass of a few 10^ M© . In TDG-N and TDG-S the 
warm gas is less than 0.1 % of the cold gas mass inferred 
from ^^CO (1-0) observations (Braine et al. 2001). The 
warm gas fraction is very small compared to the typically 
^ few percent and higher found in starbursts (Rigopolou 
et al. 2002). 

3.2. Large Scale Distribution of Non-Stellar Emission 

Figure 8 shows the stunning false color image of 
NGC 5291. The 3.6 & 4.5 /im stellar emission is coded 
blue and green, respectively, and non-stellar 8/im emis- 
sion is coded red. TDG-N and TDG-S stand out promi- 
nently in red. There is an excellent correspondence be- 
tween the 8/um sources and the 24 near-UV bright knots 
discovered by L79, with all the sources in their Figure 
1, which overlap with our areal coverage, detected. Sim- 
ilarly, all of the emission line sources detected by L79, 



Maza et al. (1991) and DM98 are visible. Figure 9 shows 
the individual IRAC band images. The non-stellar emis- 
sion in these objects almost certainly arises from strong 
7.7 /xm and 8.6 /xm PAH features that fall within the 
IRAC 8.0 fj,m filter. It is likely that any additional TDG 
candidates lie within the projected Ehi = 1 M© pc~^ 
contour overlaid on the dust image (Figure 1). The IRAC 
flux densities for the sources labelled in Figure 1 are listed 
in Table 4. A substantial fraction exist in close groups 
(e.g., #49-51), while knots #36,39,41-43 are arranged 
along a 30" (9 kpc) arc-like distribution. The bright- 
est PAH emitters are clearly associated with regions of 
higher Shi in Figure 1, though, as noted by both L79 and 
M97, there is a tendency to avoid the actual peak. Also 
note that a clump of faint knots surrounded by exten- 
sive and diffuse non-stellar emission lies 80" northwest 
of NGC 5291. A number of sources beyond the projected 
1 Mq/pc^ hi distribution also exhibit bright non-stellar 
emission. Because of the width of the 8.0 iim IRAC filter, 
knots with bright non-stellar emission but lacking an op- 
tical spectral identification that connects them with the 
system, may be at substantially different red-shifts and 
not associated with the NGC 5291 system (e.g., DM98 
measured red-shifts of 0.07 and 0.2 for sources #54 & 
#45 in Figure 1). Further observations will be required 
to confirm their physical connection with the NGC 5291 
system. 

3.3. Candidate Tidal Dwarf Galaxy Morphology 

In non-stellar emission the knots range in size from es- 
sentially unresolved (D < 2.5" or ~0.8 kpc) to -10" (3 
kpc) in extent. A number of them show interesting struc- 
ture. Sources #16 and #17, for example, appear to be 
connected by a ribbon- like bridge in Figure 10, while 
#47, 49-51 form a close group of irregularly shaped blobs. 
L79 noted that faint tails could be observed extending up 
to '--^20" from several of the brighter knots. These are 
visible in Figure 10 for TDG-N (#26) & TDG-S (#33). 
TDG-S (#33) possesses the longest tail at -38" (11 
kpc) in projection, with a noticeable brightening in its 
southernmost extent which curves off to the East. Sev- 
eral large arcs or loops can be seen in the non-stellar im- 
ages (see Figures 1, 10 & 11). Two are visible —20 and 
40" southeast of NGC 529 1's nucleus and may represent 
tidal arms emanating from this galaxy. The expected 
stellar component may be blended with the diffuse emis- 
sion. Additional wispy non-stellar emission is apparent 
~ 1.3 ' south of NGC 5291 that may be a complex tail 
from the close pair #47 & 51. While the knots invari- 
ably have lower fiux densities at 3.6 /xm and 4.5 /xm, the 
more luminous non-stellar sources are also detected in 
these bands. For example, the heads and tails in TDG- 
S (#33) and TDG-N (#26) are detected in both stellar 
3.6 /im and non-stellar 8 /im emission, as can be; scion in 
Figure 10. However, significant differences can be seen 
in the stellar and non-stellar structure of knots like #47, 
49, 50 - 51 and # 16 & 17. The tails and disturbed mor- 
phologies may be signatures of merger events between 
companion knots. 

3.4. NGC 5291 and The Seashell Galaxies 

While originally classified as an elliptical, NGC 5291 
was re-classified as an SO galaxy by L79 on the basis of 
weak spiral structure and a possible dust lane visible in 
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near-ultraviolet light. Both are consistent with our as- 
sessment that the 3.6 and 4.5 fim bands are dominated 
by cool stars, as in normal galaxies. Figure 11, top-left, 
shows a closeup of NGC 5291 and the Seashell at 3.6 fim 
using a logarithmic stretch, and with contours in Figure 
11, top-right. NGC 5291 has a bright compact stellar 
core. A slight twisting of the outer isophotes is appar- 
ent, but otherwise the stellar distribution is quite regular. 
The situation is very different in non-stellar emission, as 
shown in Figure 11, bottom- left (also with a logarithmic 
stretch), where NGC 5291 is dominated by a bright and 
marginally resolved nuclear source (R ~ 2", peak Fg.o^m 
= 225 /iJy arcsec"^) and an asymmetric 20" (6 kpc) 
diameter ring (Fg.o^tm ^ 36 /iJy arcsec"^). The ring is 
highly non-uniform, with two hot-spots north and west 
of the nucleus and a bifurcation to the southeast. The 
intensity of the non-stellar emission is also noticeably 
weaker south and west of the nucleus. Line emission de- 
tected by DM98 at a position 8" north of NGC 5291's 
nucleus coincides with this ring. No line emission was de- 
tected 4" south of the nucleus where Figure 11, bottom- 
left, shows the non-stellar emission to be considerably 
weaker. From this we deduce that the ring is actively 
forming stars, and the non-stellar emission reflects strong 
PAH features. Whether it is the result of a bar induced 
resonance (Buta & Combes 1985) or a collision as in the 
Cartwheel ring galaxy (Higdon 1996) can be determined 
in principle by detecting expansion using high angular 
resolution HI observations. 

The two faint arcs of non-stellar emission extending 
20" south of NGC 5291's nucleus in the bottom-left panel 
in Figure 11 resemble tidal arms seen in other interact- 
ing galaxies. There is in addition to these a counter- 
spray of dust projecting ^^40" (12 kpc) to the northwest 
from NGC 5291's nucleus. This can be easily seen in a 
smoothed version of the non-stellar emission (5" fwhm 
kernel) in the bottom-right panel of Figure 11, where it 
coincides with Ehi> 5 Mq pc~^ gas in Figure 1 and grad- 
ually merges into a fainter and more diffuse non-stellar 
distribution. Knots of non-stellar emission are also vis- 
ible in this material, though they are not as bright at 
8.0 /xm as the majority of the candidate TDGs. The 
extent of this material suggests a direct connection be- 
tween NGC 5291 and the TDGs and candidates even in 
the absense of the HI distribution. 

The Seashell displays a highly perturbed stellar mor- 
phology at 3.6 ^m (Figure 11, top-left), similar to that 
seen in optical light. This consists of a prominent nucleus 
and edge-on disk plus two structures extending above 
the plane, one forming an apparent ring. In non-stellar 
emission, only the nucleus and edge-on disk components 
remain, with both much fainter than NGC 5291. This is 
consistent with its lack of measured HI and line emission 
(M97; DM98). Also noteworthy is the disk galaxy #5 in 
Figure 1 (# 435 Richter 1984, M97), located ~4' south- 
west of NGC 5291. Its red color in Figure 8 indicates 
that it is rich in PAHs with a flocculent spiral structure. 
This galaxy was classified as an SA(s)c by M97. However 
a close inspection of the IRAC images shows evidence of 
a stellar bar and a dusty resonance ring. It shows no 
signs of either a recent tidal interaction or ram pressure 
stripping. 



3.5. Mid-infrared Spectral Energy Distribution 

We now have a detailed understanding of the mid- 
infrared properties of both TDG-N and TDG-S, which 
are very similar. Their IRAC SEDs are "notched" , i,e, 
the flux density falls from 3.6 to 4.5 fim and then steeply 
rises from 4.5 to 8.0 /zm. This spectral shape is char- 
acteristic of star forming regions. For an example, see 
the average 2-12 fim Infrared Space Observatory PHT- 

5 spectrum from observations of forty normal galaxies, 
shown in Figure 3 of Lu et al. (2003). As we are search- 
ing for additional candidate TDGs, we use the TDG-N 
spectrum as our template star forming galaxy. To char- 
acterize all the sources listed in Table 5, each SED is 
compared to a series of template spectra, all normalized 
at 3.6 /xm. Including TDG-N the templates are an SO 
(Devriendt, Guiderdoni & Sadat, R. 1999) an elliptical 
(Silva et al. 1998) and the narrow line Seyfert 1, I Zw 1 
from Weedman et al. (2005) as an AGN template. The 
3.6-5 /im portion of the AGN template is an extension 
of the continuum to the IRS 5-38 fim spectrum. These 
templates encompass the observed range of IRAC SEDs 
and the closest match is used to classify each source in 
Table 4. For TDG-N and TDG-S we measured the near- 
infrared J, H and Ks flux densities using a 19.5" diameter 
aperture on the 2MASS images (Cutri et al. 2003). For 
TDG-S (# 33) we quote a 5-cr upper hmit for the H-band 
magnitude. To match the extended wavelength cover- 
age of these two TDGs we included the ISO-SWS spec- 
trum of the nucleus of the starburst galaxy M 82 (Lutz 
et al. 1998) with additional J, H and Ks photometry 
from the 2MASS images matched to the ISO-SWS aper- 
ture. TDG-S shows evidence of a 1.6 fim "bump" from 
evolved stars, whereas the near infrared emission is flat- 
ter in TDG-N and a closer match to M 82, which is shown 
as a dotted line Figure 12 (dark blue line in electronic 
edition). The near- infrared component is likely from a 
combined population of evolved and low-mass main se- 
quence stars. Figure 12 shows example matches for all 
our templates^. There are 15 objects within the ring (in- 
cluding TDG-N) whose SEDs are weU fit by TDG-N (# 
15, 17, 18, 19, 25, 26/TDG-N, 32, 33/TDG-S, 36, 39, 
41, 42, 47, 49, 50). These are our candidate TDGs. An 
additional 7 sources in the ring resemble the TDG-N tem- 
plate, but with weaker PAH emission (# 5, 13, 22, 28, 
30, 43, 51). These sources are identified as star forming 
regions in Table 4. Source # 5 is the spiral galaxy SE of 
NGC 5291 coincident with the HI ring. Our results show 
that 61 % of the knots identified in the HI ring have mid- 
infrared SEDs consistent with a young stellar population 
with 41 % identified as candidate TDGs. If the spectral 
shape of the candidate TDGs match that of TDG-N and 
TDG-S both at shorter and longer wavelengths, then the 
near-infrared will probe the evolved and low-mass main 
sequence stellar population and the mid-infrared will re- 
veal a population of sources lacking a dominant cool 40 
— 60 K dust component. 

Three sources outside the ring (# 45, 54 & 59) are a 
good match to the TDG-N template. These include the 
two known background starburst galaxies in DM97 (# 54 

6 #45). Two sources (#64, 66) are a reasonable match 
to TDG-N. Being outside of the projected HI, these five 

^ SED matches for the full sample are given at 
|http://isc. a stro . Cornell . edu/"sjuh/tdg/NGC5291/ 
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sources are listed as star forming regions. Some sources 
are a poor match to TDG-N but have "notched" SEDs 
indicating some star formation. NGC 5291 's (#44) SED 
matches the SO template, while the Seashell's (#35) SED 
is consistent with both an elliptical and SO template. 
Source #20 is an elliptical galaxy. A number of sources, 
#14, 23, 24, 38, 52, 53, 57, 60, 65, 69, may be elliptical 
galaxies or foreground stars, and #2, 9, 10, 21, 61 may 
be SO galaxies or foreground stars. These are labelled E* 
and S* respectively, in Table 4. Outside of the projected 
HI ring there are a 7 sources that may be background 
AGN as their SEDs can be fit with a simple power law 
(#3, 4, 6, 8, 12 56 73). 

The number of sources located inside and outside the 
projected HI ring are roughly equal, but the ring has only 
~ 40 % the areal coverage of the background. 44 % of 
the sources located within the ring arc good matches to 
the TDG-N SED compared to 8 % of the sources located 
outside the ring. To estimate the sample contamination 
we assume that the three good matches to TDG-N that 
lie outside the HI ring are from foreground/background 
sources. From this number we infer that approximately 
two of the candidate TDGs may not be associated with 
the NGC 5291 system. 

3.6. Color Magnitude and Color-Color Relations 

In this section wc look at the color-color and color- 
magnitude plots to see whether TDG-N, TDG-S and 
the thirteen candidate TDGs occupy a distinct region 
of color-color/magnitude space when compared to other 
galaxy types. In Figure 13 we plot the [4.5] — [8.0] ver- 
sus [3.6] — [4.5] color. The vertical and horizontal lines 
mark the color zero-points, using the zero-point flux den- 
sities listed in §§ 2.2. The plot is now divided into four 
quadrants. A source with a flat SED would lie at the 
crossing point. Sources in the lower left quadrant have 
falling SEDs typical of E and SO galaxies like the SeashcU 
and NGC 5291 (sec # 35 & 44 in Figure 12). AGN-hke 
sources, for example, #12a shown in Figure 12, with ris- 
ing SEDs are located in the upper right quadrant. The 
lower right quadrant contains sources with active star 
formation. The candidate TDGs, for example, see #26, 
33 & 42 in Figure 12, are plotted with a framed-diamond 
symbol. They have [4.5] - [8.0] > 3.2 and [3.6] - [4.5] 
< 0.4. Sources which are a reasonable match to TDG- 
N, for example see #13 in Figure 12, but with weaker 
8 /im emission and therefore not included in our candi- 
date TDG list, are plotted as double diamonds. Sources 
which are a good match to TDG-N, but lie outside the 
projected HI ring, are plotted as solid diamond symbols. 
The remaining "notched" SED sources in this quadrant, 
for example, see #48a in Figure 12, are a poor match to 
TDG-N and are shown as open-diamond symbols. 

For comparison the sample of 18 galaxies from Pahre et 
al. (2004) are plotted as crosses. This sample contains 
representative galaxies from the Hubble morphological 
sequence and spans El through IB(s)m galaxies. The 
candidate TDGs do not overlap with this sample. The 
late type galaxy NGC 5669 has the reddest [4.5] - [8.0] 
color of 2.7 in the Pahre et al. (2004) sample, as com- 
pared to the average value of 3.8 for the candidate TDGs. 
Also included in Figure 13 are the 19 BCD galaxies from 
the KPNO International Spectroscopic Survey (KISS), 
observed with IRAC by Rosenberg et al. (2005). These 



sources, shown as filled circles, were selected from the 
KPNO International Spectroscopic Survey (KISS, Salzer 
et al. 2000) conducted in the NDWFS Bootes field (Jan- 
nuzi & Dey 1999), and span a wide range in metallicity 
(7.8 < 12 -f logO/H < 9.1) and star formation rate (0.1 < 
SFR < 2.1 M© yr^i). Only one BCD (#2316) overlaps 
with the TDG and candidate TDG distribution. This 
object is noteworthy for having the highest oxygen abun- 
dance and second highest SFR (1.3 Mq yr~^) among the 
BCDs. The two next closest BCDs in Figure 13 (#2328 
& #2359 in Rosenberg et al. 2005) have the highest and 
third highest SFRs. Otherwise, their properties (i.e.. Ha 
equivalent widths, metallicity) do not distinguish them 
from the remaining BCD sample. The TDGs and can- 
didate TDGs all show significantly enhanced non-stellar 
emission, most likely due to PAHs, relative to normal 
spirals and even BCD galaxies. 

In Figure 14 we plot the [3.6] — [8.0] color as a func- 
tion of the absolute [8.0] magnitude. There is no overlap 
with the Rosenberg et al. (2004) sample of BCDs. The 
average for the two TDGs and the 13 candidate TDGs is 
^1.5 magnitudes redder in [3.6] — [8.0] color and ^1.9 
magnitudes fainter at 8 /im. The TDGs and candidate 
TDGs have a very narrow range in [3.6] — [8.0] color of 
0.7 compared to the BCD sample, which has a spread of 
2.5. 

3.7. Stellar Masses For The Candidate TDGs 

A widely used technique for determining the stellar 

mass distributions in disk galaxies is to multiply the mea- 
sured broadband luminosity by the appropriate mass- 
to-luminosity ratio (van Albada et al. 1985). This is 
generally best carried out at near-infrared wavelengths, 
where the eff'ects of dust and massive stars are mini- 
mized (Kranz et al. 2001; Verhcijen 2001). Adopting 
this approach, we estimated the stellar masses for TDG- 
N & TDG-S by first calculating their luminosities in the 
Kg-band, defined as Lk^ = Aird'^uF^ with v=lA x 10^"' 
Hz. Expressed in terms of the sun's Kg-band luminosity 
(1/0, Ks = 5.3 X 10^^ W), this can be written as 

Lk, = 3.3 X 10-37 d^ Fk3 (Lq, kJ, (7) 

where d is the distance in meters and Fk^ is the fiux 
density in Janskys. For TDG-N & TDG-S we found 
= (3.4 and 6.3) x 10^ L©, jf^, respectively. Stellar 
masses were obtained by multiplying these luminosities 
by Kg-band mass-to-luminosity ratios (M^/Lxg) calcu- 
lated by Bell & de Jong (2001) for model spiral galaxies. 
This choice seemed appropriate given that both the metal 
abundances (Z ^ 1/3 Tjq) and optical colors (B-R ^ 0.6) 
of the candidate TDGs in the NGC 5291 system are simi- 
lar to those found in the disks of late spiral galaxies (L79; 
DM98). Bell & de Jong found that their derived M*/L 
depended on the optical color, as was expected. However, 
the color dependence is minimized in the near-infrared. 
For a B - R = 0.6, their models give M*/Lk = 0.3-0.6 
M0/L0, K^- Using these values results in a stellar mass 
of 1.5 + 0.6 X 10^ M© for TDG-N, which is consistent 
with Bournaud et al. 's (2004) estimate obtained from its 
blue luminosity and a M/Lb of two. For TDG-S we de- 
rive a stellar mass of 3.0 1 1.0 x 10*^ Mq. These are 
roughly 10% the stellar mass of the LMC (Leroy et al. 
2005). 
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We also estimated stellar masses for the thirteen can- 
didate TDGs in Table 4 showing evidence of strong PAH 
emission like TDG-N & TDG-S. These are the framed- 
diamond symbol sources in Figures 13 and 14. Since 
none of these sources were detected by 2MASS at J, H, 
or Ks, we calculated Fk, flux densities from their 3.6 
/xm data using the mean (Kg — [3.6]) color index for 
TDG-N & TDG-S, which is 0.9 t 0.5. For these sources, 
Fk^ in the above equation was replaced with Fs.e^m x 
j^0-o.4xo.9 X (645./277.5), where the last factor takes 
into account the different zero points in the two bands. 
We derive a mean Kg luminosity and stellar mass of 4 
t 2 X 10"^ L0, and 2 t 1 X 10^ Mq respectively for 
these sources, with the largest source of uncertainty aris- 
ing from the observed spread in Kg — [3.6] color. TDG 
candidate #32 ('^l' southeast of TDG-N) is the the most 
massive of this group with M, = 6 1 3 x 10^ M©. De- 
fined this way, TDG-N & TDG-S possess stellar masses 
that are roughly an order of magnitude larger than the 
average for the thirteen other candidate TDGs of the HI 
ring with [4.5]-[8.0] > 3.2. 

In a similar way, we estimated stellar masses for the 
nineteen BCD galaxies from Rosenberg et al. (2005). 
Four of the BCDs were detected at Kg as well, and their 
average Kg — [3.6] color (1.0 0.5) is similar to the mean 
for TDG-N & S, with a comparable spread in value. For 
these galaxies we either calculated Lq^ and di- 
rectly from their Kg-band flux densities, or indirectly us- 
ing their 3.6 /xm flux densities and mean Kg — [3.6] color. 
For the latter fifteen objects, the largest source of imcer- 
tainty in the Kg — [3.6] colors, which we take to be the 
observed spread among the four BCDs detected at Kg- 
band. On average, the BCDs have considerably larger 
stellar masses than the TDG-candidate population as- 
sociated with NGC 5291, with a median M* = 6^3 
X 10^ M0. TDG-N & TDG-S are close to this average 
value, however the other TDG-candidates have stellar 
masses roughly an order of magnitude smaller than the 
typical BCD in Rosenberg et al.'s KISS sample. These 
objects arc detected to considerably larger distances than 
NGC 5291, leading to a possible Malmquist bias. How- 
ever, we find a nearly identical median stellar mass for 
the Rosenberg et al. sample if we restrict ourselves to 
those BCDs within 125 Mpc, which is the distance where 
our faintest TDG-candidate would have been detected in 
all four IRAC bands at 5 a. Of the sources in Table 4, 
only TDG-N & TDG-S appear to have stellar masses ex- 
pected for a typical BCD. 

4. CONCLUSIONS 

We have presented mid-infrared spectra of two tidal 
dwarf galaxies, TDG-N & TDG-S in the NGC 5291 sys- 
tem. Both are PAH rich with features at 6.2, 7.7, 8.6, 
11.3 and 16.5 /xm. The 5-14 /xm spectra are remark- 
ably similar in shape to both the starburst galaxy M 82 
and the reflection nebula NGC 7023. However, detailed 
examination reveals their differences. The relative PAH 
strengths are indicative of an equal mix of ionized and 
neutral PAHs, with groups of ~ 100 carbon atoms, as op- 
posed to the almost totally ionized PAHs present in the 
starburst galaxy M 82. The rise in flux density (F^) flat- 
tens off around 20 /xm in both TDG-N and TDG-S and 
the spectra can be fit with a T ~ 140 K dilute blackbody. 



This is in marked contrast with many spiral and starburst 
galaxies, such as M 82, whose spectra peak around 60 - 
100 /um due to a dominant cooler dust component with a 
temperature ^ 40 K. In TDG-N the 1.6 fim bump from 
an evolved stellar population is absent. Both TDG-N 
and TDG-S are dominated by young stellar populations. 
Their spectra exhibit many fine structure and molecular 
lines allowing us to calculate the physical properties of 
the star forming regions. The electron density derived 
from the [Slll] ratio in TDG-N is 600 < Ue < 3000 
cm-3. Both the [Slv]/[Siii] and [Ne iii]/[Ne ii] ra- 
tios are higher than typical ratios measured in starburst 
galaxies and closer to the values measiirecl for BCD and 
Wolf-Rayet galaxies. Both metallicity and the age of a 
starburst effect the level of excitation. The TDGs' have 
one-third solar metallicity and the [Ne iii]/[Ne ii] line ra- 
tios are consistent with models that have the most recent 
epoch of star formation occurring less than five million 
years ago. This results in a moderately high excitation. 
The age of the starburst is consistent with models which 
predict that the young stars have formed in-situ and were 
not part of the parent galaxy. 

For the first time we have the sensitivity to measure 
the low J rotational lines from molecular hydrogen and 
thereby make a direct measurement of the warm molec- 
ular mass. Assuming the gas is in LTE with an ortho- 
to-para ratio of 3, we derive temperatures of 460 ^ 75 K 
and 400 1 40 K, with a corresponding mass of warm gas 
of (1.0 t 0.7) and (1.5 t 0.6) xlO^ M© in TDG-N and 
TDG-S, respectively. This warm mass is < 0.1 % of the 
cold gas mass inferred from ^^CO (1-0) measurements. 

IRAC observations have identified 13 additional 
sources within the HI ring that resemble TDG-N and 
TDG-S. These candidate TDGs occupy a distinct region 
of IRAC color space with [3.6] - [4.5]' < 0.4 and [4.5] - 
[8.0] > 3.2 as compared to elliptical, spiral, irregular and 
BCD galaxies, which tend to have [4.5] - [8.0] < 3.2. The 
combined TDG and candidate TDG sample has a very 
narrow dispersion in [3.6] — [8.0] color of ~ 0.7, unlike a 
sample of BCDs, which has a spread of ~ 2.5. Compared 
to the BCDs the average candidate TDG is ~ 1.5 magni- 
tudes redder with an [3.6] — [8.0] color of 3.6 and ~ 1.9 
magnitudes fainter at 8 fim, with an absolute magnitude 
of ~ -21.6. Approximately 40 % of the knots identified 
in the HI ring fit the TDG template with strong PAH 7.7 
& 8 fxm emission a factor ~ 6 times the luminosity of the 
cool stellar population. Using their 2MASS Kg-band flux 
densities and a starburst M^/Lk^ appropriate for their 
B-R colors, we derive stellar masses of (1.5 1 0.6) and 
(3.0 t 1.0) X 10^ M0 for TDG-N and TDG-S respec- 
tively. These values are comparable to similarly defined 
stellar masses for a sample of BCDs, though still only 
~10% of the LMC's stellar mass. Adopting a mean Kg 
— [3.6] color of 0.9 for the remaining candidate TDGs, 
we used their 3.6 /im luminosities to derive stellar masses 
from 0.8-6.0 x 10^ Mq. 

The IRAC images of the HI ring show an incredibly 
rich nursery for star formation. Many of the candidate 
TDGs have a disturbed morphology possibly resulting 
from previous tidal encounters with fellow companions. 
These are new star forming systems and may be local 
analogues of star formation resulting from tidal interac- 
tions in the early universe. 
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Future work will compare this population of TDGs 
with our Spitzer sample of TDG systems, which includes 
observations of Stephan's Quintet, Arp 245 and Arp 105. 
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Foundation; the SIMBAD database; and archival data 
from the National Radio Astronomy Observatory, which 
is a facility of the National Science Foundation operated 
under cooperative agreement by Associated Universities, 
Inc. 
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TABLE 1 
Observational Parameters 



Object RA (J2000) Dec (J2000) Date Instrument Integration Execution Time 





h:m:s 


d:m:s 






(s) 


(min) 


NGC5291 


13:47:23.00 


-30:25:30.00 


2004-02-17 


IRAC 


3 X 12 


44.87 


TDG-N 


13:47:20.50 


-30:20:51.00 


2004-07-17 


IRS 
IRS-SLl 
IRS-SL2 
IRS-SH 
IRS-LH 


60 X 4 
60 X 4 
60 X 4 
120 X 7 


77.81 


TDG-S 


13:47:23.00 


-30:27:30.00 


2004-07-17 


IRS 
IRS-SLl 
IRS-SL2 
IRS-SH 
IRS-LH 


60 X 4 
60 X 4 
60 X 4 
120 X 7 


77.73 



Spitzer TDGs 



TABLE 2 
PAH Flux 



Keutrai Flux EW PAH^/PAH7.7 



(^tm) 


/in— 21 ■\x7„vvi — 2"^ 
(lU Wcm ) 






■T^ l~\ ' AT 








a o 

D.Z 


1 Q O ~l~ 1 A 


-o.y _ u. * 


U.OO _ U.UD 


'7 '7 


'i'7 o "I" r; o 
o/.z _ O.Z 


-Iz.U 1.4 


i.UU 


8.6 


8.1 1.3 


-2.1 0.4 


U.zz _ U.Uo 


11.3 


8.lj^ 0.6 


-1.7 ^ 0.2 


U.Zz _ U.Uo 


12.0 


1.8l|] 0.6 


-0.4^ 0.2 


A nc; "T n ao 
U.Uo _ U.Uz 


12.6 


5.7 ^ 1.2 


-1.3 0.3 


0.15 _ 0.04 


13.5 


1.5 ^ 0.5 


-0.4 ^ 0.1 


0.U4 _ U.Ol 


16.5 ^ 


0.6 ^ 0.2 


-U.Ul _ U.Ul 


O.Oz _ 0.01 


TDG-S 








6.2 


3.3 t 1.0 


6.8 1 2.1 


0.33 1 0.11 


7.7 


10.0 t 1.3 


12.3 1 1.8 


1.00 


8.6 


2.6 i 0.7 


2.5 t 0.7 


0.27 1 0.08 


11.3 


2.3 t 0.1 


1.54 t 0.30 


0.23 t 0.05 


12.0 


0.6 t 0.2 


0.4 t 0.2 


0.06 t 0.03 


12.6 


1.2 t 0.4 


0.8 1 0.3 


0.12 t 0.04 


13.5 


0.10 t 0.06 


0.1 1 0.1 


0.01 1 0.01 


16.5 ^ 


1.0 t 0.5 


-0.03 t 0.02 


0.10 t 0.05 



''Prom the IRS-SH aperture, uncertainty includes aperture scaling 
factor between IRS-SL and IRS-SH. 



TABLE 3 
Line Fluxes 





TDG-N 




TDG-S 




Line 


Line Flux 


EW 


Line Flux 


EW 




(10 Wcm-2) 


(^tm ) 


(10 Wcm-2) 


(^jm ) 


[S IV] 10.51 Mm 


7.40 1 0.15 


-O.OlSl 0.002 


1.05 1 0.22 


-O.OO2I 0.001 


[Ne II] 12.81 lira 


10.65 t 0.06 


-0.0241 0.001 


4.05 1 0.31 


-O.Olll 0.003 


[Ne III] 15.56 lira 


25.31 1 1.03 


-0.065t 0.002 


5.47 1 0.11 


-0.0161 0.002 


[S III] 18.71 /im 


15.45 t 0.58 


-0.04lt 0.002 


4.66 1 0.17 


-0.0151 0.000 


[Siii] 


10.61 t 4.02 


-O.lOOl 0.049 


<37.87 t 0.54 




0-0 8(3) 9.66=^ 


2.2 t 0.9 


-0.04 1 0.02 


1.6l0.5 


-0.08 1 0.03 


0-0 8(2) 12.28 


1.9 1 0.8 


-0.0031 0.001 


0.9 1 0.4'= 


-O.OO2I0.OOI 


0-0 S(l) 17.03 


1.1 1 0.4 


-0.0031 0.001 


1.3 1 0.3<i 


-O.OO2I 0.001 


0-0 S(0) 28.22 


< 1.7^ 




<1.7« 





^IRS-SL 

''FWHM fixed to instrumental resolution of 0.088 /im. 
'^FWHM fixed to instrumental resolution of 0.020 lira. 
''FWHM fixed to instrumental resolution of 0.028 //m. 
''3-C7 upper limit. Calculated using a 0-order fit to the continuum. 
The upper Umit is 3 x RMS x FWHM. The FWHM is 0.05 //m. 
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TABLE 4 

IRAC Flux Densities in an 8 Pixel Diameter Aperture 



Object^- RA (J2000) Dec (J2000) Fa.g^™ F4 









(°) 


(°) 




(i"Jy) 


(i"Jy) 




1 




AGN 


206.77913 


-30.463076 


82.613.5 


67.613.4 


73.2I15.2 


67.ll 10.3 


2b 




SO* 


206.78200 


-30.411720 


220.0^3.5 


2OI.0I3.4 


104.0I15.2 


IO7.0I1O.3 


3" 




AGN 


206.78366 


-30.425823 


114.0l3.5 


I23.0I3.4 


167.0ll5.2 


I47.0I1O.2 


4b 




AGN 


206.80237 


-30.392924 


126.0l2.6 


I75.0I2.9 


268.01 10.9 


365.0l9.2 


5cd 




SF 


206.80253 


-30.456886 


5460.0j:20.2 


339O.0I18.I 


6200.01 65.4 


16500.01 55.3 






AGN 


206.80267 


-30.439257 


188.0^3.5 


343.0I3.4 


6I4.0I15.2 


1240.01 10.3 


7b 




N 


206.80351 


-30.363409 


52.0^3.5 


33.2I3.4 


7O.1I14.9 


62.31 10.3 


8b 




AGN 


206.80364 


-30.372922 


66.813.5 


97.1I3.4 


2O6.0I15.2 


240.01 10.3 


9 




SO* 


206.80908 


-30.468743 


IO7.0I3.5 


74.8I3.4 


79.51 15.2 


61.61 10.3 


IQb 




SO* 


206.81047 


-30.420336 


169.0l3.5 


I22.0I3.4 


9O.8I15.2 


77.0I1O.3 


11 




AGN? 


206.82215 


-30.339189 


112.0l3.5 


II2.0I3.4 


3OI.0I14.3 


I6I.0I1O.O 


12b 




AGN 


206.82216 


-30.414722 


35.2I2.6 


44.0I2.9 


6O.2I1O.9 


IOO.0I9.2 


13 


H - 


SF 


206.82294 


-30.472132 


140.0^3.5 


IO7.0I3.4 


I7O.0I15.2 


419.0I1O.2 


14 




E* 


206.82355 


-30.453369 


204.0I3.5 


I27.0I3.4 


7I.6I15.2 


51. ol 10.2 


15 


- f 


TDG 


206.82406 


-30.366257 


39.913.5 


29.1I3.4 


86.61 15.2 


252.01 10.2 


16 




N 


206.82418 


-30.356122 


145.0I3.5 


I22.0I3.4 


92.31 15.2 


23I.0I1O.2 


17 


- h 


TDG 


206.82499 


-30.352538 


89.913.5 


57.1I3.4 


259.0I15.2 


663.01 10.2 


18 




TDG 


206.82504 


-30.361186 


45.0I3.5 


24.813.4 


I3O.0I15.2 


352.0I1O.3 


19 




TDG 


206.82672 


-30.344119 


30.713.5 


9.613.4 


< 36.3 


I88.0I1O.2 


20b 


1 C 


E 


206.82682 


-30.383323 


1050. 0l6. 5 


625.0l6.2 


34I.0I28.O 


246.01 18.5 


21 




SO* 


206.82883 


-30.452879 


207.0I3.5 


I33.0I3.4 


93.71 15.2 


75.0I1O.2 


22 


- b 


SF 


206.83111 


-30.448559 


49.5I2.6 


32.1I2.9 


40.71 10.8 


IO6.0I9.2 


23 




E* 


206.83137 


-30.415181 


243.0I3.5 


167.0l3.4 


II8.0I15.2 


73.0I1O.3 


24 




E* 


206.83159 


-30.399068 


47I.0I3.5 


3I6.0I3.4 


224.0I15.2 


101. ol 10.2 


25 


Aj 


TDG 


206.83236 


-30.339936 


69.3I2.6 


38.0I2.9 


402.0I1O.8 


486.0l9.2 


26<i<= 


B i 


TDG 


206.83508 


-30.347620 


513.0I6.5 


389.316.2 


1844.5128.2 


4804.01 18.8 


27 




N 


206.83516 


-30.447661 


8O.7I2.6 


56.612.9 


28.91 10.8 


76.419.2 


28 




SF 


206.83619 


-30.445187 


I9.8I2.6 


8.1I2.9 


18.91 10.9 


7I.1I9.3 


29 






206.83759 


-30.384396 


<11.7 


<9.0 


26.91 10.8 


85.2I9.2 


30 


4 - 


SF 


206.83873 


-30.472597 


373.0I3.5 


256.0l3.4 


286.01 15.2 


I35O.0I1O.3 


31 






206.84210 


-30.383212 


8.9I2.6 


IO.4I2.9 


52.31 10.9 


158.0l9.2 


32 


- g 


TDG 


206.84312 


-30.362525 


II8.0I3.5 


72.5I3.4 


3OO.0I15.2 


853.01 10.2 


33<^f 


F a 


TDG 


206.81521 


-30.458553 


295.0I6.5 


160.916.1 


617.8128.2 


1930.01 18.8 


34b 




AGN? 


206.84561 


-30.491370 


92.7I3.5 


8O.8I3.4 


78.21 15.2 


111. ol 10.3 


35sb 


- B 


E SO 


206.84707 


-30.417140 


I75OO.0I2I.2 


10400.01 19.9 


8050.01 72.0 


6IIO.0I54.7 


36 




TDG 


206.84741 


-30.449865 


27.4I2.6 


17.8I2.9 


28.41 10.9 


174.0I9.2 


37b 




AGN? 


206.84765 


-30.489807 


6I.1I2.6 


59.6I2.9 


48.71 10.8 


72.1I9.2 


38 




E* 


206.84767 


-30.357765 


895.0l3.5 


578.0I3.4 


432.0I15.2 


227.0I10.3 


39 




TDG 


206.84773 


-30.447670 


I7.1I2.6 


I5.0I2.9 


27.31 10.8 


II3.0I9.3 


40b 




AGN? 


206.84823 


-30.345555 


5I.4I3.5 


79.7I3.4 


52.31 15.2 


101. ol 10.2 


41 


- c 


TDG 


206.84868 


-30.445188 


25.0I2.6 


IO.8I2.9 


4O.0I1O.9 


I37.0I9.2 


42 


E - 


TDG 


206.84913 


-30.452056 


5I.0I3.5 


32.7I3.4 


IO5.0I15.2 


258.01 10.2 


43 


D d 


SF 


206.84915 


-30.439878 


62.8I2.6 


4O.5I2.9 


62.1I10.9 


I8O.0I9.2 


44gi 


- A 


E 


206.85203 


-30.407169 


69400.0l41.9 


40800.01 36.1 


30900.0l87.2 


36400.0l88.7 


45b 


- 2 


SF 


206.85413 


-30.374597 


76.0I2.6 


62.0I2.9 


25.1I1O.8 


198.0l9.2 


46b 






206.85697 


-30.492522 


282.0l3.5 


176.0l3.4 


71.91 15.2 


60.91 10.3 


47 




TDG 


206.85752 


-30.430211 


38.612.6 


17.3I2.9 


97.61 10.8 


284.01 9.2 


48b 




N 


206.85858 


-30.395310 


205.0I2.6 


I8I.0I2.9 


149.0I1O.9 


29O.0I9.2 


49 




TDG 


206.85934 


-30.433200 


27.5I2.6 


I5.1I2.9 


5O.1I1O.8 


158.0l9.2 


50 




TDG 


206.86013 


-30.434832 


27.7I2.6 


2I.0I2.9 


57.2I10.8 


I62.0I9.3 


51 


C e 


SF 


206.86030 


-30.430790 


I34.0I2.6 


88.9I2.9 


142.01 10.8 


313.0I9.2 


52b 




E* 


206.86159 


-30.508598 


215.0I3.5 


175.0I3.4 


I2I.0I15.2 


76.71 10.3 


53b 




E* 


206.86201 


-30.359309 


268.0I3.5 


173.0I3.4 


IO4.0I15.2 


75.2I1O.3 


54b 


2 1 


54 


206.86375 


-30.371721 


663.0l3.5 


446.0l3.4 


629.0I15.2 


4070.01 10.3 


55b 






206.86543 


-30.447166 


II5.0I3.5 


III.0I3.I 


IO4.0I15.2 


64.91 10.2 


56b 




AGN 


206.86557 


-30.418478 


78.612.6 


88.4I2.9 


83.41 10.9 


2O8.0I9.2 


57b 




E* 


206.86574 


-30.367851 


173.0I3.5 


II6.0I3.4 


73.0I15.2 


58.61 10.2 


58b 




AGN? 


206.86612 


-30.414799 


14.6I2.6 


I9.4I2.9 


17.81 10.8 


15O.0I9.2 


59b 




SF 


206.87124 


-30.476743 


236.0l3.5 


156.0l3.4 


252.0I15.2 


944.0I1O.2 



Spitzer TDGs 



TABLE 4 

- CONT. 





Object^ 




RA (J2000) 
(°) 


Dec (J2000) 
(°) 


(MJy) 


^ 4.5fim 


(i"Jy) 




F8.0(jm 

(MJy) 




60'' 


3 - 


E* 


206.87193 


-30.404070 


809.0I3.5 


480.0l3.4 


293.0I15, 


.2 


222.0I10, 


.3 


ei** 




SO* 


206.87499 


-30.370020 


407.0^3.5 


306.0^3.4 


258.0^15, 


.2 


221. olio. 


.3 






N 


206.87787 


-30.404201 


95.5+2.6 


86.7+2.9 


42.1+10, 


.9 


124.0+9, 


,2 


63*' 




N 


206.87844 


-30.455269 


80.0"'"2.6 


68.9+2.9 


<16, 


.8 


162.0+9, 


.3 


64^ 




SF 


206.88901 


-30.390435 


133.0+3.5 


117.0^3.4 


96.0+15, 


.2 


284.0+10, 


.3 


65^ 




E* 


206.89271 


-30.364745 


923.0I3.5 


650.0^3.4 


404.0ll5, 


.2 


261. olio. 


.3 


66*" 




SF 


206.89412 


-30.467691 


69. 3^3. 2 


64. 5^3. 3 


31.1^14 


.9 


135. olio. 


,2 


67*= 






206.89467 


-30.431946 


119.0t3.3 


82.9^3.2 


87.91 15, 


.2 


111. olio. 


,2 


68'' 




N 


206.89631 


-30.437985 


6O.2I2.6 


53.1I2.9 


28.51 10, 


.9 


87.7l8, 


.8 


69'' 




E* 


206.89914 


-30.395302 


201.0^3.5 


133.0t3.3 


IO2.0I15, 


.2 


54.61 10. 


.3 


70^ 






206.90090 


-30.367889 


241.0^3.5 


158.0^3.3 


I24.0I15, 


.2 


175.0I1O, 


.2 


71^ 




N 


206.90105 


-30.460684 


243.0^3.5 


151.0^3.4 


75.6I15, 


.2 


I3I.0I1O, 


.2 


72^ 






206.90456 


-30.421585 


37.8j:3.5 


49.913.3 


36.51 14, 


.4 


II4.0I1O, 


,2 


73.^ 




AGN 


206.91035 


-30.448771 


45.7^3.5 


43.5t3.4 


82.61 15, 


.2 


94.71 10, 


,2 


74b 






206.91079 


-30.450288 


22.7^3.2 


65.0^3.3 


<41, 


.2 


73.5I9, 


.7 


75^ 




N 


206.91286 


-30.439546 


142.0^3.5 


I29.0I3.4 


36.41 15, 


.2 


204.0I1O, 


.0 



^Second column is the L79/M97 and DM98 naming convention. 
Third column is the classification from the SED match in SS3.5, T; 
candidate TDG, SF: reasonable fit to TDG-N, star-forming region, 
N: "notched" SED but poor match to TDG-N, E: elliptical, E*: 
elliptical/foreground star, SO : SO, SO* : SO/foreground star, AGN 
:AGN, AGN? : not a simple power-law 

''Source is outside HI ring. 

<^Galaxy 435 (Richter 1984, M97). 

''I6 pixel diameter aperture. 

«TDG-N. Source K09.09 (Pena, Ruiz & Maza 1991) 
'^TDG-S., the tail is source G in L79 

^Mag. auto flux density as source is extended beyond 16 pixel 
diam. aperture. Note the aperture correction for an infinitely 
extended source has not been applied. The correction factors listed 
in the SOM are 0.94, 0.94, 0.63, 0.69 at [3.6] , [4.5] , [5.8] , [8.0] , 
respectively. 

•'Seashell. 

'NGC 5291. 
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Fig. 1. — The IRAC "non-stellar" image of the NGC 5291 system 
created by subtracting stellcir light from the 8.0 //m image (see 
text). This emission is most likely dominated by strong 7.7 fj,m 
and 8.6 jum PAH bands. Arrows identify sources listed in Table 
4. For example, # 26 & 33 are TDG-N and TDG-S, respectively. 
NGC 5291 is # 44 and the Seashell is #35. The contours show 
the distribution of HI (Shi = 1, 3, 5, 7 & 9 Mq pc"^) created by 
reprocessing the VLA DnC archive data (Sfwhm = 50"). Original 
VLA data published by Malphrus et al. (1997). 
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Fig. 2.— Spitzer IRS spectra of TDG-N. (Top) IRS-SL spectrum 
showing strong PAH features and 0-0 8(3) 9.7 /Lim molecular 
Hydrogen. (Bottom). Combined IRS-SH & IRS-LH spectrum 
showing PAHs and fine structure lines. The flattening of the 
continuum emission beyond A > 20 ixm signifies the absence of 
a dominant cool 40 K dust component. 
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Fig. 3.— Spitzcr IRS spectra of TDG-S. (Top) IRS-SL spectrum 
showing strong PAH features and 0-0 S(3) 9.7 /im molecular 
Hydrogen. (Bottom). Combined IRS-SH & IRS-LH spectrum 
showing PAHs and fine structure lines. The flattening of the 
continuum emission beyond A > 20 fj,m. signifles the absence of 
a dominant cool 40 K dust component. 
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Fig. 4. — A comparison between the IRS Short-Low spectrum of 
TDG-N (thick black Une), the ISO-SWS spectrum of the starburst 
galaxy M 82 dash line (dashed red-line electronic edition) from 
Lutz et al. (1998) after smoothing to the IRS Short-Low resolution 
and normalized to 7 fim, and an IRS Short-Low spectrum of the 
reflection nebula NGC 7023 normalized to 7 fim thin line (dashed 
blue-line electronic edition) from Werner et al. 2004b, showing the 
remarkable similarity in shape. Note the [Ne ll] 12.81 fira line is 
blended with the 12.6 fim PAH. 
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Fig. 5. — Spitzcr spectra of emission lines detected in TDG-N 
(histogram). The hnes were fit with a Gaussian profile combined 
with a first order polynomial fit to the continuum (solid line) . The 
lines are from IRS-SH/LH, except 0-0 S(3), which was observed 
with IRS-SL. 
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Fig. 6. — Spitzer IRS spectra of emission lines detected in TDG- 
S (histogram). The Hnes were fit with a Gaussian profile combined 
with a first order polynomial fit to the continuum (solid line) . The 
lines are from IRS-SH/LH, except 0-0 S(3), which was observed 
with IRS-SL. All lines were detected in both nod positions except 
for the 0-0 S (2) line. 
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Fig. 7.— Excitation Diagrams for TDG-N (left) and TDG-S 
(right). The absolute value of the reciprocal of the slope is the 
excitation temperature. 
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Fig. 8. — IRAC false color image of the NGC 5291 system, defined 
by 3.6 /im (blue), 4.5 /im (green), 5.8 /xm (orange) and 8.0 /fm (red). 
Sources dominated by starlight (e.g, NGC 5291 &; 'Seashell') are 
blue in color while objects with strong PAH emission, including 
TDGs and the flocculent spiral galaxy (#5 in Table 4), are red in 
color. 
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Fig. 9.— (Top-left), IRAC 3.6 image of the NGC 5291 
system. A linear stretch between -0.04 and 1.6 MJy sr^-"^ is used. 
The scalebar represents 1', i.e., a linear scale of 18 kpc. (Top-right), 
IRAC 4.5 //m image. A linear stretch between -0.05 and 1.2 MJy 
sr~^ is used. (Bottom-left), IRAC 5.8 ^m image. A linear stretch 
between -0.13 and 0.71 MJy sr'^ is used. (Bottom-right), IRAC 
8.0 /im image. A linear stretch between -0.10 and 0.72 MJy sr~^ 
is used. 
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Fig. 9 — cont. 
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Fig. 10.— Close-ups of TDG-N (# 26), TDG-S (# 33) and 
candidate TDGs associated with NGC 5291 showing emission at 
3.6 //m (top-row) and "non-stellar" emission from PAHs (bottom 
row). The vertical bar represents 30", i.e., a lineaj: scale of 9 kpc. 
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Fig. 11.— IRAC images of NGC 5291 and the "Seashell" 
galaxy, showing (top-left) 3.6 /im emission, (top-right) 3.6 ;um 
emission with surface brightness contours, (bottom-left) "non- 
stellar" emission, and (bottom-right) "non-stellar" emission after 
smoothing with a 5" Gaussian kernel. The images are displayed 
with a logarithmic stretch between 0.1 and 1.0 MJy sr~i, and the 
contours are equally spaced in log-space between 0.3 and 30.0 MJy 
sr~l. 
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Fig. 12. — Example IRAC Spectral Energy Distributions (SEDs) 
overlaid with a series of template galaxies: TDG-N the solid line 
(colored pink in electronic edition), Elliptical the dot-dash line 
(Silva et al. 1998, colored green in electronic edition), SO the dot- 
dot-dot-dash line (Devriendt, Guiderdoni & Sadat, R. 1999, colored 
red in electronic edition), M 82 nuclear region the dotted line (Lutz 
et al 1998, colored dark blue in electronic edition) and an I Zw 1 
the dash-line (Weedman et al. 2005, colored light blue in electronic 
edition). Sources with the suffix 'a' are outside the projected HI 
ring. TDG-N is #26, TDG-S is #33. #42 is an example of a good 
match to TDG-N i.e., a candidate TDG, #13 is a reasonable match 
to TDG-N and characterized as a star forming region, #48a is 
a "notched" SED indicating star formation, the Seashell (#35) 
matches both the elliptical and SO templates, NGC 5291 (#44) 
matches the SO template and #12a is a power-law. characteristic 
of an AGN. 
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Fig. 13.— [4.5] -[8.0] versus [3.6] - [4.5] color - color plot. 
The horizontal and vertical lines mark the color zero-points. A 
source with a flat SED would lie on the crossing point of these 
two lines. Sources in the lower left quadrant have falling SEDs 
typical of NGC 5291 and the Seashell, i.e., E and SO galaxies. 
AGN like sources, with rising SEDs, are located in the upper right 
quadrant. The lower right quadrant contains sources with active 
star formation. TDG-N and TDG-S are shown as a small solid star 
and square, respectively. The diamonds represent the remaining 
sources from Table 4. In particular, framed-diamonds denote the 
13 candidate TDGs. Sources located within the HI ring that are 
a reasonable match to TDG-N, but with weaker PAH emission, 
are shown as double-diamonds. Sources outside the HI ring with 
TDG-like SEDs are shown as solid diamonds. Source identification 
numbers are from Table 4. For comparison the Pahre et al. (2004) 
galaxy sample are shown as crosses and the Rosenberg et al. (2005) 
BCDs are shown as solid circles. The TDGs and candidate TDGs 
all show significantly enhanced non-stellar emission, most likely 
due to PAHs, relative to normal spirals and even BCD galaxies. 
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Fig. 14. — [3.6] - [8.0] color versus [8.0] absolute magnitude. 
Sec the caption to Figure 13 for symbol definitions. The average 
positions of the TDG/TDG-candidate and BCD samples are 
indicated by the large star and circle, respectively. The candidate 
TDGs occupy a distinct region, being on average ~ 2.2 magnitudes 
redder in [3.6] — [8.0] color and 2.2 magnitudes fainter at 8 /im than 
the Rosenberg et al. (2005) BCD sample. 
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